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The hadronic k ± -spectrum and the gluon to quark average multiplicity ratio r — N g /N q inside a high en- 
ergy sub-jet are determined from a precise definition of the jet axis, including corrections of relative magnitude 
O(y / o7) with respect to the Modified Leading Logarithmic Approximation (MLLA), in the limiting spectrum ap- 
proximation (assuming an infrared cut-off Qo — ^- QCD )• The results for the fc^-spectrum in the limiting spectrum 
approximation are found to be, after normalization, in impressive agreement with measurements by the CDF 
collaboration and the ratio has not been measured yet. 



1. Introduction 

The collimation of hadrons in jets is a ba- 
sic phenomenon of high energy collisions and its 
quantitative understanding is an important task 
for QCD. Simple differential characteristics of a 
jet are their energy and multiplicity angular pro- 
files. The collimation of energy and multiplicity 
in the jet follows from the dominance of gluon 
bremsstrahlung processes in the parton cascade 
evolution. Whereas the former is more sensitive 
to the hard processes inside a jet the second one 
is sensitive to the soft parton emissions from the 
primary parton. 

The characteristics of soft particle production, 
such as particle multiplicities, inclusive distri- 
butions and correlation functions, are derived 
in QCD in the Double Logarithmic Approxima- 
tion (DLA) and the Modified Leading Logarith- 
mic Approximation (MLLA) (for review, see pQ), 
which takes into account the leading double log- 
arithmic terms and the single logarithmic cor- 
rections. These azimuthally averaged quantities 
can be obtained from an evolution equation for 
the generating functional of the parton cascade. 
This equation provides also Next-to-MLLA (NM- 
LLA) corrections taking into account energy con- 
servation of parton splittings with increased ac- 
curacy. The corresponding hadronic observables 
can be obtained using the concept of Local Parton 



Hadron Duality (LPHD) [2] , which has turned out 
a successful description of many hadronic phe- 
nomena (see, e.g. [3]). 

The inclusive k ± -distribution of particles inside 
a sub-jet with respect to the jet axis has been 
computed at MLLA accuracy in the limiting spec- 
trum approximation [4], i.e. assuming an infrared 
cutoff Qo equal to A qcd (A = ln(Q /A qcD ) = 0) 
(for a review, see also [T]). MLLA corrections, 
of relative magnitude O (^/a~s) , were shown to be 
quite substantial [I]. Therefore, corrections of 
order 0(a s ), that is NMLLA, have been incor- 
porated in [516] . 

Another characteristic prediction of QCD is the 
asymptotic increase of the mean average particle 
multiplicity in a gluon jet over the quark jet by 
the ratio of color factors r — N g /N q — > N c /Cf = 
9/4, which is also modified by the similar set of 
corrections 0(^/a^) [71819) . In particular, this 
analysis has been extended to multiparticle pro- 
duction inside sub-jets with a precise prescription 
of the jet axis, which has been, for both observ- 
ables, identified with the direction of the energy 
flux [TU]. 

The starting point of this analysis is the MLLA 
Master Equation for the Generating Functional 
(GF) Z = Z(u) of QCD jets [1], where u = u{k) 
is a certain probing function and fc, the four- 
momentum of the outgoing parton. Together 
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with the initial condition at threshold, the GF de- 
termines the jet properties at all energies. Within 
this logic, the leading (DLA, 0(^fb~2)) and next- 
to-leading (MLLA, 0(a s )) approximations are 
complete. The next terms (NMLLA, 0(a^ 2 )) are 
not complete but they include an important con- 
tribution which takes into account energy conser- 
vation and an improved behavior near threshold. 
Indeed, some results for such NMLLA terms have 
been studied previously for global observables and 
have been found to better account for recoil ef- 
fects and to drastically affect multi-particle pro- 
duction |ll|12j . 

The main results of this work have been pub- 
lished in [516110] . Experimentally, the CDF 
collaboration at the Tevatron reported on k ± - 
distributions of unidentified charged hadrons 
in jets produced in pp collisions at ^Js = 
1.96 TeV [13]. 

2. MLLA evolution equations 

2.1. Inclusive spectrum 

We start by writing the MLLA evolu- 
tion equations for the fragmentation function 
Dg (x/z, zEOq, Qo) °f a parton B (energy zE 
and transverse momentum k ± — zEQq) into a 
gluon (represented by a hadron h (energy xE) 
according to LPHD [5j) inside a jet A of energy 
E for the process depicted in Fig. [T] As a conse- 




A = (Q, G) 



Figure 1. Parton A with energy E splits into par- 
ton B (respectively C) with energy zE (respec- 
tively, (1 — z)E) which fragments into a hadron h 
with energy xE. 



quence of angular ordering in parton cascading, 
partonic distributions inside a quark and gluon 
jet, Q,G(z) = x/z D'^ G (x/z,zEQ ,Qo), obey 
the system of two coupled equations [14] (the sub- 
script y denotes d/dy) following from the MLLA 
master equation [T] 



Qy = J da^Sf (*) [(Q(l -z)-Q)+ G{z) 
^ g (z)(G(z) + G(l-z) 



G 



dz- 



G 



(z)(2Q(z) 



G 



(1) 



(2) 



$^(z) denotes the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) splitting functions [1 , 



T R [z 1 + (1 - zY 



&Jz) = 2N c (- + cf£(z) 



(3) 



where 4> 9 q (z) = z-2 and 4> 9 g (z) = (z- 1)(2- z{l - 
z)) are regular as z -> 0, C F = {N% - 1)/2N C , 
Tr = H//2 (N c — 3 is the number of colors for 
SU(3) C and n/ = 3 is the number of light flavors 
we consider). The running coupling of QCD (a s ) 
is given by 

2ir 11 4 

a s = ; r, 4N C 3 = —N c - ~T R 

and 

£={l/x) , y = ln(k x /Q ) , A = ln(Q /A QCD ), 

where Qo is the collinear cut-off parameter. 
Moreover, 



G = G(l) = xD^(x,E@o,Q ), 
Q = Q(l) = xD^(x,EQ ,Qo). 



(4) 



At small x <C z, the fragmentation functions be- 
have as 



B ( z ) ~ Pb In-, In— — 

x Qo 



P h B (]nz + £,y), 



p'g being a slowly varying function of two log- 
arithmic variables hi(z/x) and y that describes 
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the "hump-backed" plateau [15]. Since recoil ef- 
fects should be largest in hard parton splittings, 
the strategy followed in this work is to perform 
Taylor expansions (first advocated for in [16]) of 
the non-singular parts of the integrands in (|1I2[) 
in powers of lnz and ln(l — z), both considered 
small with respect to £ in the hard splitting region 
z ~ 1 - z = O (1) 



B(z) = B{\)+B t (l)\nz + (in 2 z) 



(•5) 



with z 1 — z. Each ^-derivative giving an ex- 
tra v / o7 factor (see [E]), the terms Bg(l) In z and 
1^(1) In (1 - z) yield NMLLA corrections to the 
solutions of {5}. Making used of ([5]), after in- 
tegrating over the regular parts of the DGLAP 
splitting functions, while keeping the singular 
terms unchanged, one gets after some algebra 
( 7o 2 = 2N c a s /Tr) [1 



Q(£,y) = S(£) + ^ [ d£' fdy'^{£' + y') x (6) 

1 - ~ ai 5{l' -£) + a 2 S(£' - t)^t(£' , v')] G(jt, y'), 

G(l,y) = 5(1) + f d£' [ y dy'^(e' + y')x 
Jo Jo 

1 - ai S(£' - £)a 2 S(£' - £)Mt\ v')] G(t, y'), (7) 

where the constants take the values at sw 0.935, 
oi = 3/4, a 2 « 0.08 and ~a 2 ~ 0.42 for n f = 
3. Defining F(£,y) = y%(£ + y)G(£,y), we can 
exactly solve the self-contained equation by 
performing the inverse Mellin transform: 



F(£,y) 



dudv 



e ui e vy F(u,v), 



(2wi) 2 

such that the NMLLA solution reads 



(8) 



G(£,y) = (£ + y + X) 



s) 



dtodv 
(2m) 2 



ds 



(lu + s) V 



where 



and 



CT 



f3 (uj - v)' 



0"1+CT 2 



Ol 

Po 



(9) 
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As can be seen, the NMLLA coefficient a 2 is 
very small. Therefore, the NMLLA solution © 
of |3 can be approximated by the MLLA solu- 
tion of G(£,y) (i.e. taking a 2 ~ 0), which is 
used in the following to compute the inclusive k ± - 
distribution. As demonstrated in [T3] , taking the 
limits a 2 sa and A w in ©, the integral rep- 
resentation can be reduced to the known MLLA 
formula, which can be written in terms of Bessel 
series in the limiting spectrum approximation [1] . 
To get a quantitative idea on the difference be- 
tween MLLA and NMLLA gluon inclusive spec- 
trum, the reader is reported to the appendix B of 
[5]. The magnitude of a 2 , however, indicates that 
the NMLLA corrections to the inclusive quark jet 
spectrum may not be negligible and should be 
taken into account. After solving ([7]), the solu- 
tion of ([6]) reads 

Q(i,y) = [G(£,y) + (ai - a^G^y) 

+ (at (ai - a^j + a 2 - a 2 ) G u (£, y) (10) 

It differs from the MLLA expression given in [3] 
by the term proportional to Gu, which can be 
deduced from the subtraction of (Cf/-/V c )x to 
Eq. ©. 

2.2. Particle mean average multiplicity 

Integrating the system (|H2p over the energy 
fraction x leads the average multiplicity of par- 
ticles produced in the jet: 



N a (Yq 



dxDA(x, E®Qi Qo)- 



where Y@ = £ + y — ln(EO /Q ). Therefore, an 
equivalent system of two-coupled evolution equa- 
tions can be deduced and expanded, following 
the same logic that led to (|6I7[) . However, in 
this case, the expansion in y/a s is performed for 
Yq 3> lnz ~ ln(l — z), where similarly, z ~ 1. We 
limit ourselves here to the NMLLA expression of 
r that reads 



— = rn(l 



njo 



r 2ll) 



(11) 



where the asymptotic value of r is r = N c /Cf = 
9/4. The MLLA term r\ has been calculated in 
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[7]. The coefficients Tk can be obtained from the 
Taylor expansions of N\ (Yq + Inx) and N\ (Y& + 
ln(l - a;)) for Y e > lnx and Y e > ln(l - x) in 
the evolution equations [16] . The values of for 
n f =3 read r x = 0.185 and r 2 = 0.426 16J. The 
NMLLA solution for the mean multiplicity in a 
gluon jet is found as [TB] 



NZ(Y e )~)C(Ye) 



exp 



-Cl/00 




2c 2 



YZ 



(12) 



e 



with ci = 0.28, C2 = 0.38 for nj = 3 and /C, the 
LPHD normalization factor. The pre-exponential 
term (Ye)~ Cl//3 ° is the MLLA contribution to N£, 
while the one oc C2, the NMLLA one. 



Ft 



( duD% Q (u,EQ ,uEQ) 

l J X 

xD^(^-,uEe,Qo), (14) 



A=g,q 



u being the energy fraction of the intermediate 
parton A. E\ q describes the probability to emit 
A with energy uE off the parton Aq (which ini- 
tiates the jet), taking into account the evolution 
of the jet between 6o and <d. D\ describes the 
probability to produce the hadron h off A with 
energy fraction x/u and transverse momentum 
k ± ps uEQ > Qo (see Fig. [2]). As discussed in 



3. Single inclusive k± distribution of 
charged hadrons and the gluon to quark 
mean average multiplicity ratio r at 
NMLLA 

3.1. Inclusive k± distribution 

Computing the single inclusive kj_— distribu- 
tion and the ratio r = N g /N q inside a sub-jet re- 
quires the definition of the jet axis. The starting 
point of our approach consists in considering the 
correlation between two particles (hi) and (h2) of 
energies E\ and E^ which form a relative angle 
inside one jet of total opening angle 6 > & [TT] . 
Weighting over the energy E 2 of particle (h2), 
this relation leads to the correlation between the 
particle (h=hl) and the energy flux, which we 
identify with the jet axis (see Fig. [2]) [4]. Thus, 
the correlation and the relative transverse mo- 
mentum k±_ between (hi) and (h2) are replaced 
by the correlation, and transverse momentum of 
(hi) with respect to the direction of the energy 
flux. Finally, we obtain the double differential 
spectrum d 2 A r /dxd6 of a hadron produced with 
energy E\ = xE at angle (or k± w xEQ) with 
respect to the jet axis. As demonstrated in [4], 
the correlation reads 



d 2 A^ 



dzdlne dln6 



FJJ (x,e,E,Q ), 



(13) 



where F'i 

Ao 

lution of 



is given by the convo- 
two fragmentation functions 




(Jet Axis) 



Figure 2. Inclusive production of hadron h at 
angle O inside a high energy jet of total opening 
angle 0o and energy E. 



[4], the convolution (fl"4|) is dominated by u ~ 1 
and therefore D^ o (u, EOq, uEQ) is determined 
by the DGLAP evolution pQ. On the contrary, 
the distribution D\ (^,uEQ,Q ) at low x <C u 
reduces to the hump-backed plateau, 

D\ (l + In u, y) X ^ U p\ (I + In u, Y @ + hxu), (15) 

with Y = £ + y = ]n(EQ/Q Q ). Performing the 
Taylor expansion of D to the second order in 
(In u) and plugging it into Eq. (fl4|) leads to 



Fg, q {^y) = jf{C) g ,S,y) G(l,y). 



(16) 



where (C) g ^ q can be written in the symbolic way 
{C) m {l, y) ~ f {l, y)+fi(£, y)V^ s +f 2 (e, y)a„, (17) 
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where fk follow from the DGLAP evolution equa- 
tions pQ and can be symbolically written as 

/ fc ~ f duu(\nu) k Diju,EQ a ,uEQ) (18) 

J X 

and corrections C(a™^ 2 ) from the n— logarithmic 
derivatives of the inclusive spectrum: 

D A 8£ n - W{as >■ 

The exact formulae of the color currents are re- 
ported in |4|5|6j . Indeed, since soft particles are 
less sensible to the energy balance, the correla- 
tion (|14[) disappears for these particles, leading to 
the sequence of factorized terms written in (I17p . 
The first two terms in Eq. (|17p correspond to the 
MLLA distribution calculated in [4] when D\ is 
evaluated at NLO and its derivative at LO. NM- 
LLA corrections arise from their respective calcu- 
lation at NNLO and NLO, and, mainly in prac- 
tice, from the third term, which was exactly com- 
puted in |5)6j . Indeed, since x/u is small, the 
inclusive spectrum D\(£,y) with A = G, Q are 
given by the solutions |9]) and (fT0|) of the next- 
to-MLLA evolution equations ^ and (O respec- 
tively. However, because of the smallness of the 
coefficient G(£,y) shows no significant differ- 
ence from MLLA to NMLLA. As a consequence, 
we use the MLLA expression (|9]) for G(£, y) with 
a 2 = 0, and the NMLLA (TTUl) for Qi^y)- The 
functions Fg and are related to the gluon dis- 
tribution via the color currents (C) g . q which can 
be seen as the average color charge carried by the 
parton A due to the DGLAP evolution from Aq 
to A. 

This calculation has also been extended beyond 
the limiting spectrum, A =/= 0, to take into ac- 
count hadronization effects in the production of 
"massive" hadrons, m = O (Qo) [18]. We used, 
accordingly, the more general MLLA solution of 
([7]) with ai = for an arbitrary A ^ 0, which 
can only be performed numerically. The NM- 
LLA (normalized) corrections to the MLLA re- 
sult are displayed in Fig. [3] for different values 
A = 0,0.5, 1. It clearly indicates that the larger 
A, the smaller the NMLLA corrections. In par- 
ticular, they can be as large as 30% at the lim- 
iting spectrum (A = 0) but no more than 10% 




Figure 3. NMLLA corrections to the color current 
of a quark jet with Yq =6.4 and £ = 2 for various 
values of A. 



for A = 0.5. This is not surprising since A ^ 
(Qo ^ A qcd ) reduces the parton emission in 
the infrared sector and, thus, higher-order cor- 
rections. 

The double differential spectrum (d 2 N/d£dy) 
can now be determined from the NMLLA color 
currents (| 17[1 by using the quark and gluon dis- 
tributions. Integrating over £ leads to the sin- 
gle inclusive y-distribution (or k ± -distribution) of 
hadrons inside a quark or a gluon jet: 



dN 
Ay 



Ye -v 



d( 



d 2 N 
d£dy 



(19) 



A 



The MLLA framework does not specify down to 
which values of I (up to which values of x) the 
double differential spectrum (d 2 A/d£dy) should 
be integrated over. Since [<d 2 N jdldy) becomes 
negative (non-physical) at small values of I (see 
e.g. [4]), we chose the lower bound £ m ; n so as to 
guarantee the positiveness of (d 2 N jd£dy) over 
the whole £ m - m < £ < Yq range (in practice, 
£ 9 min ~ 1 and £^ in ~ 2). Having successfully com- 
puted the single k ± -spectra including NMLLA 
corrections, we now compare the result with exist- 
ing data. The CDF collaboration at the Tevatron 
recently reported on preliminary measurements 
over a wide range of jet hardness, Q — EQq, 
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Figure 4. CDF preliminary results for the inclu- 
sive k ± distribution at various hardness Q in com- 
parison to MLLA and NMLLA predictions at the 
limiting spectrum (Qo = Aqcd); the boxes are 
the systematic errors. 



in pp collisions at */s = 1.96 TeV [13]. CDF 
data, including systematic errors, are plotted in 
Fig. 0] together with the MLLA predictions of 
[I] and the present NMLLA calculations, both 
at the limiting spectrum (A = 0) and taking 
A = 250 MeV. The experimental distribu- 

QCD r 

tions suffering from large normalization errors, 
data and theory are normalized to the same bin, 
ln(fcj_/l GeV) = —0.1. The agreement between 
the CDF results and the NMLLA distributions 
over the whole k ± -range is particularly good. In 
contrast, the MLLA predictions prove reliable in 
a much smaller k ± interval. At fixed jet hard- 
ness (and thus le ), NMLLA calculations prove 
accordingly trustable in a much larger x interval. 

Finally, the kt-distribution is determined with 
respect to the jet energy flow (which includes 
a summation over secondary hadrons in energy- 
energy correlations). In experiments, instead, the 



jet axis is determined exclusively from all parti- 
cles inside the jet. The question of the matching 
of these two definitions goes beyond the scope of 
this letter. The NMLLA k ± -spectrum has also 
been calculated beyond the limiting spectrum, 
by plugging ((9]) with a 2 = into (fT9|) . as illus- 
trated in Fig. El However, the best description of 
CDF preliminary data is reached at the limiting 
spectrum, or at least for small values of A < 0.5, 
which is not too surprising since these inclusive 
measurements mostly involve pions. Identifying 
produced hadrons would offer the interesting pos- 
sibility to check a dependence of the shape of k ± - 
distributions on the hadron species, such as the 
one predicted in Fig. [5] Moreover, the softening 
of the k x -spectra with increasing hadron masses 
predicted in Fig. [5] is an observable worth to be 
measured, as this would provide an additional 
and independent check of the LPHD hypothesis 
beyond the limiting spectrum. This could only 
be achieved if the various species of hadrons in- 
side a jet can be identified experimentally. For- 
tunately, it is likely to be the case at the LHC, 
where the ALICE [TJ] and CMS [5D] experiments 
at the Large Hadron Collider have good identi- 
fication capabilities at not too large transverse 
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momenta. 

3.2. Ratio r = N g /N q 

Integrating (jT5J) over the energy fraction x 
yields the corresponding sub-jet multiplicity 
of hadrons inside the angular range < 60 of 
the jet A [10] 

N% o (e;E,G )n I duuDi Q (u,EQ ,uEQ) 

A=q,g Ju « 

xN%(uEe,Q ), (20) 

where N% is the number of hadrons (partons at 
scale Qo) produced inside the sub-jet A of total 
virtuality uEQ > Q and u — Q /EQ. Within 
the leading parton approximation, we expand the 
multiplicity N^(uEQ,Q ) in {20]) at u ~ 1. Sim- 
ilarly to the logic applied in (|15p for the compu- 
tation of the k±— distribution, we take the log- 
arithmic dependence of the average multiplicity 
N\, such that for In it < Y® = \n(E@/Q ) and 
EQ ^> A, one obtains the average multiplicity 
ATjj (Ye , Yq) of soft hadrons within the sub-jet 
of opening angle with respect to the energy flow 

N^ g (Y eo ,Y e ) « ^(C) q , g (Y @0 , Y e )N^(Y e ), (21) 

where N^Yq) is written in (fT2"]) and 
(C) qi g(YQ , Yq) is the average color current of 
partons forming the energy flux. We write the 
color current as 

(C) g ,g(Ye ,Ye) ~ f + fiVa s + f 2 a s , 

where ff. = fk (Y<~> , Yq ) follow from (JTTJ) and 
(US]) and corrections 0(a^ 2 ) from the logarith- 
mic derivatives of the average multiplicity (|12[) : 

N A dYg ~ s ■ 

The exact formulae of (C) g)g (Ye , lg>) are given 
in the appendix of [10 . The ratio of the gluon to 
the quark jet average multiplicity reads 




(C) g (Y eo ,Y e ) 
(C) q (Y eo ,Y e y 



For this quantity, in the limit — > ©o, the ap- 
propriate ratio (JTTJ) r = jf- = r (l - ri7 - r 2 7o) 
is recovered. 



Next we study the consequences for the behav- 
ior of the ratio at full angles ©o ~ 1 and for sub- 
jets at reduced angles 0. Results for the ratio 
using our formulae for the color currents (C) q ^ g 
[TO] are depicted in Fig. [6] The full line shows the 
ratio (fTTj) for an isolated Ao jet as a function of 
Yq — \n(EQ /Qo)- The three curves below cor- 





NMLLA 




E=50GeV 




E=100GeV 




E=1000GeV 







3 4 5 6 7 8 



Y 8 



Figure 6. Multiplicity ratio N g /N q for jets and 
sub-jets as function of energy variable Y = 
ln(£6/Qo). 



respond to the sub-jet energies E = 50, 100, 1000 
GeV and the intersection points with ([TT| to the 
limit — ► 0o- Therefore, the mixing of quark 
and gluon jets reduces the multiplicity in the 
gluon jet and increases it in the quark jet at suf- 
ficiently high energies and the effect from inter- 
mediate processes in the new definition of multi- 
plicity ratio amounts to about 20% at a reduced 
energy scale Yq ~ Ye — 2. 

4. Conclusion 

We have studied a definition of the jet which 
arises from the 2-particle correlation. The jet axis 
corresponds to the energy weighted direction of 
particles in a given cone ©o- In this way smaller 
angles (or transverse momenta) can be meaning- 
fully determined. 
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The new effect from 2-particle correlations is 
the appearance of the "color current" in the re- 
sult, which reflects the possibility of intermediate 
quark- gluon processes. We then studied the con- 
sequences of this analysis on the measurement of 
differential (fc^-distribution) and global (particle 
multiplicity) observables, where we have added 
MLLA and NMLLA corrections to the known LO 
results pp. 

The single inclusive k ± -spectra is determined 
within this logic and the agreement between NM- 
LLA predictions and CDF data in pp collisions at 
the Tevatron [13] is impressive, indicating very 
small overall non-perturbative corrections and 
giving further support to LPHD [2] . 

The main for the ratio r = N g /N q reflects the 
possibility of intermediate quark-gluon processes. 
This effect can be seen by the scale breaking be- 
tween jets and sub-jets at the same scale EQ but 
different energies E and opening angles O. Typi- 
cal effects are of the order of 20%. However, these 
effects on this observable have not been measured 
yet. 
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